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Abstract 

Corrections to scaling, associated with deviations of the order parameter 
from the scaling morphology in the initial state, are studied for systems with 
0(n) symmetry at zero temperature in phase-ordering kinetics. Including 
corrections to scaling, the equal time pair correlation function has the form 
C(r, t) = fo{r/L) + L^^ fi{r / L) + • • •, where L is the coarsening length scale. 
The correction-to-scaling exponent uj and the correction-to-scaling function 
fi{x) are calculated for both nonconserved and conserved order parameter 
systems using the approximate Gaussian closure theory of Mazenko. In gen- 
eral u; is a non-trivial exponent which depends on both the dimensionality, 
d, of the system and the number of components, n, of the order parameter. 
Corrections to scaling are also calculated for the nonconserved 1-d XY model, 
where an exact solution is possible. 



Typeset using REVTgX 



I. INTRODUCTION 



The dynamics of a system undergoing phase-ordering following a quench from the high 
temperature (disordered) phase to the ordered phase is of great interest [0]. The kinet- 
ics of systems with 0{n) symmetry subject to 'Model A' dynamics [0(i.e. systems with 
nonconserved order parameter) and 'Model B' dynamics 0] (systems with conserved order 
parameter) have been previously studied within a Gaussian closure theory originally de- 
veloped by Mazenko |^,^ following the seminal work of Ohta, Jasnow and Kawasaki (OJK) 
1^. In previous work |Q we have computed the form of the corrections to the scaling limit, 
and the correction-to-scaling exponent, for a number of systems with nonconserved order 
parameter. These include some exactly soluble models, and the Model A dynamics of a 
scalar field within the Mazenko theory. 

In the present work we turn our attention to systems with a vector order parameter, 
both nonconserved and conserved. The corrections to scaling for systems with continuous 
symmetry will be calculated using the Mazenko theory. It should be mentioned that this 
approach has been shown to be more successful, at a quantitative level, in systems with 
nonconserved order parameter than those with conserved order parameter . Nevertheless, 
the results obtained in the conserved case are in qualitative agreement with those obtained 
in simulations. Furthermore, the Mazenko approach seems the only available method to 
probe the questions of corrections to scaling addressed here. In particular, we found that, 
for nonconserved scalar fields, the correction-to-scaling exponent u is predicted by this 
approach to have a nontrivial value. We will show that this same feature is present for the 
vector fields, with and without conservation. 

It is well established [0 that at late times most phase-ordering systems approach a scaling 
regime, where the equal-time pair correlation function C{r, t) = (0(x+r, t) •0(x, t)) takes the 
form C{r,t) = f[r/L{t)]. The characteristic length scale L{t) grows with time as L{t) ~ t", 
where a is the growth exponent which depends on the nature of the dynamics and the 
symmetry of the order parameter. In particular, a = 1/2 for nonconserved order parameter 
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systems, while a = 1/3 for systems with conserved scalar order parameter and a = 1/4 
for systems with conserved vector order parameter (with logarithmic corrections for n = 2, 
d > 2 [|1^). In previous work ^ we studied how scaling is approached in nonconserved 
order parameter models such as the 1-d Ising model with Glauber dynamics, the n— vector 
model with n = oo, the approximate OJK theory and the Mazenko theory for scalar fields. 
In all these cases u was found to be trivial {u = 4) except the last, for which u was found 
to be non-trivial and dimensionality dependent. The relevance of corrections to scaling lies 
in interpreting experimental and simulation results, where it is advantageous to know how 
the scaling limit is approached. Corrections to scaling in systems with finite n > 1 in d = 3 
and d = 2 were not considered in The main objective in this article is to study systems 
with n > 2. 

This article is devoted to the study of the corrections to scaling for systems with 0{n) 
symmetry in phase-ordering dynamics. The leading corrections to scaling enter the correla- 
tion function in the form 

C{r,t) = fo{r/L) + L-Mr/L), (1) 

where fo{x) is the 'scaling function' and fi{x) the 'correction-to-scaling function'. The 
quantity which unites theory, computer simulation and experiment is the structure factor 
S{k,t) = L'^gQ^y) + L'^~'^ gi{y) , where go{y) and gi{y) are the d-dimensional Fourier trans- 



forms of fo{x) and fi{x) respectively, and y = kL. Coniglio and Zannetti |[TT| solved the 
conserved 0{n) model for n = oo exactly, and found that no simple scaling exists. Instead 
a 'multiscaling' behavior was obtained, raising the question of whether simple scaling exists 
in conserved order parameter systems with n > 1 generally (or even for conserved scalar 
fields). However, it was later shown by Bray and Humayun [|1^], analytically within the 
Mazenko theory, that scaling does exist for large but finite n. Attempts to find multiscaling 
behaviour in simulation date for conserved scalar fields [|13|, or the conserved XY model 
in two |[T^ or three [|T^] dimensions were not successful. It is now generally believed that 



scaling is recovered asymptotically in time in the conserved 0{n) model, for all finite n, 



though multiscahng may be observable in the preasymptotic regime 



There are other sources of corrections to scahng apart from the one considered in this 
paper. In phase-ordering systems there is, in addition to the time-dependent coarsening 
scale L{t), a second characteristic length scale - the 'defect core size' ^ - in systems with 
topological defects. The corrections to scaling associated with nonzero defect core size (where 
^ is the domain wall thickness in scalar systems) are expected to enter as a power of C,/L. 
Here we are interested primarily in the corrections to scaling associated with non-scaling 
initial conditions. We therefore suppress the contributions associated with nonzero core size 
^ by taking the 'hard-spin' limit, i.e. working with an order-parameter field whose length is 
everywhere unity, 0^ = 1, which forces ^ = (though in the Mazenko theory this limit will be 
taken at the end). Also thermal fluctuations at T > may give rise to significant corrections 
to scaling for systems quenched to a nonzero final temperature T (where < T < T^, 
with Tc the critical temperature) as has been shown explicitly in the nonconserved 0{n) 
model with n ^ oo [0. However, we will only be studying systems quenched to T = 0. 
Although corrections to scaling due to thermal fluctuations and nonzero ^ are important we 
will not consider them further in this paper. 

The outline of the paper is as follows. In the following section the approximate Mazenko 
theory is discussed and some general concepts are introduced. Section III deals with noncon- 
served order parameter systems. In section IV, corrections to scaling for the nonconserved 
1-d XY model will be studied. Systems with conserved order parameter are considered in 
section V. Section VI concludes with a summary and discussion. 



II. MAZENKO THEORY 



A 'Gaussian closure' theory, building on the earlier work of Ohta, Jasnow and Kawasaki 
has been developed by Mazenko [Q. This theory has been successfully applied to 0{n) 
models in the theory of phase-ordering dynamics . The equation of motion for an order 
parameter (p with continuous symmetry, for systems quenched to T = 0, is 
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where p = 1 and p = for conserved order parameter (Model B) and nonconserved order 
parameter (Model A) systems respectively. In (^, V"(0) is a symmetric double- well potential 
for the scalar case, and a 'wine bottle' potential with a degenerate continuum manifold for a 
vector order parameter. Compact notation has been used in which '1' represents the space- 
time point (xi,ti) and means the Laplacian with respect to Xi. Multiplying (Q) by 0(2), 
averaging over initial conditions, and using the translational invariance of C(12) gives (for 

tl=t2 = t) 
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where now is the Laplacian with respect to r = |xi — X2I and C(12) = (0(1) • 0(2)). 
The angular brackets denote the average over the initial conditions. In order to evaluate 
the average of the last term in (^) one introduces an auxiliary field m(r, t) related to by 
V^0 = 2 dV{(f))/d(j), with boundary condition m/\m\ as |m| — > 00, and = at m = 0. 
Near a defect, the field m(r) is the position vector of the point r in the plane normal to the 
defect. The assumption that m is a Gaussian field enables the evaluation of the average of 
the last term on the right hand side of (^ giving 
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25-0(1) ^7 

where Sq = (m(l)^) and 7(12) = (m(l)m(2))/[(m(l)^)(m(2)^)]^/^ is the normalised correla- 
tor of the field m (where m is one of the components of m). An explicit expression which 



relates 7 to C(12) was given in ||T9 
C 
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where B{y, z) = T{y)T{z)/T{y + z) is the Beta function and F(a, 6; c; z) the hypergeometric 
function. Equations (^ and (|^) provide closed form equations for C(12). On substituting 
(S) in (H) one obtains an equation for 7 which can in principle be solved numerically and 



substituted back into to obtain the correlation function C(12). We note at this point 
that in deriving the correlation function (^), the 'hard-spin' limit (f) = m/\m\ was employed. 
Since this result holds far from defect cores, it will correctly describe the scaling limit where 
the defects are dilute. Here we are also using it to compute the corrections to scaling. 



III. NONCONSERVED 0(N) MODEL 

For a nonconserved system p = 0, and equation (^) is simply 

For n = 1, using the properties of the hypergeo metric function the last term on the right 
hand side of @ can be written in terms of C(12) only, resulting in an equation which is 
independent of 7(12). Corrections to scaling in this case where obtained in our previous 
work [Q, and will not be considered further here. For general n, 7 cannot be eliminated 
in favour of C(12), and we will therefore work with 7 instead of C(12). From dimensional 
considerations we see that Sq ~ and can be chosen as Sq = L^/A. This choice effectively 
defines L, up to an overall constant. For n 00, an expansion in 1/n can be performed 
on C(7), and in this limit 'jdC/d'j = C + C^/n + 0(l/?7,^). For n = 00, Mazenko theory 
reduces to the n = 00 n-vector model for which an exact solution, including the corrections 
to scaling, is known |^. Expressing (|^) in terms of 7 explicitly leads to 

2 dt \dr J dr^ r dr 2L^ ^ ' ^ ^ 

where C-y = dC/d'~^ etc. Since C{r,t) is a function of 7(r, t), the scaling and corrections 
to scaling can be imposed on 7(r, t). In the scaling limit we expect 7(r, t) to approach the 
scaling function 7o(r/L) which is L-independent if all lengths are scaled by L. In this limit 
therefore one expects LdL/dt = constant. Including corrections to scaling in 7(r, t) and L{t) 
as usual ^ we can write 

7(r,t) =70 (^) +1^-^71 (^)+---, (8) 
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7=70 



and 6 is a constant. Equating leading and next-to- leading powers of L in the usual way gives 
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with = [dC/dj] 7=70 etc. The primes indicate derivatives with respect to the scaling 
variable x = r/L. 

Equations (|TB|) and ([T^) are to be integrated numerically subject to appropriate 'initial' 
conditions imposed at x = 0. Since x = corresponds to 70 = 1, the initial conditions 
are obtained by considering the regime 70 ^ 1. Using the properties of the hypergeometric 



functions one can derive relations between C(7o) and its derivatives as 70 — >■ 1. Up to 



prefactors of order unity, we find in this limit 
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and so on. We have given explicit expressions for C^f^-y^/Cyg and (^^07070/^70 as 70 — 1 for 
the values of n which we are going to study. Using the above results one can show that 
the small-x behavior of 7o(x) is given by 

7o(x) = l- Aa;2^... (Ig) 

for n > 2, where the limiting forms in (|15D were used to demonstrate that the term involving 
C70 70 7^*70 in dH) is sub dominant for n > 2. 

For large-x, 70 (also C(12) 0) and equation ( |13|) becomes linear because in 
this limit the second term in ([l3|) is negligible. It is easy to show that two linearly in- 
dependent solutions of the linearised equation have the asymptotic forms 701 ~ x'"^^ and 
7o2 ~ x^'^"'^ exp(— for x 00. As equation ([T3|) is integrated forward from x = 0, 
the large-x solution obtained will in general be a linear combination of 701 and 702- The 
amplitudes of 701 and 702, however, depend on A. For systems with initial conditions con- 
taining only short-range spatial correlations (as is the case for systems quenched from high 
temperature), a power-law decay is unphysical, and A is determined by the condition that 
the coefficient of the power-law term, 701, must vanish Note that A is related to the 
exponent A describing the decay of the autocorrelation function via X = d — X. Values 
for A are given in Table 1 for 2 < ri < 5 and 1 < d < 3. Comparison of the predicted values 
of A with simulations ||2^ and experiments p3| show reasonable agreement. It can be shown 



that for d 00 the OJK result (A = d/2) is recovered for both scalar and vector 
cases. The same limit for A is also obtained for n — ^ 00 at arbitrary d. 

The correction-to-scaling exponent, u, is found from (|1^) in a similar way to the de- 



termination of A from (0). In order to specify initial conditions for the numerical inte- 
gration of (0), we need the small-x behavior of 71 (x). A small-x analysis of (|I4|) gives 
7i = bXx'^/16d{d -|- 2) -|- ■ ■ -, where the results in (pr5| ) were used to show that the last two 
terms in ( p!^ are subdominant as x — > 0. The required initial conditions are therefore 
7i(0) = 7i(0) = 0- As X ^ 00, the last two terms in (|T3|) can be neglected. The two linearly 
independent solutions of the simplified equations have a power-law tail (~ x"*^'^'^^^^) and a 
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Gaussian tail (~ x'^ exp(— a;^/8)) for large x, where q = 2X — d + uj ii u > 2 and q = 2X — d + 2 
otherwise. Having already found A, u; is chosen on physical grounds in the same way as A, 
namely that the coefficient of the power-law term in the large-x solution should vanish. The 
values of uj obtained are given in Table 2 for 2 < n < 5 and 1 < d < 3. Note that u ^ 4 for 
d oo, as the OJK result (and its generalization to vector fields) is recovered in this limit. 

After solving (|^) and (0) for 70 (x) and 71 (x) we use these results to get the scaling 
function /o(x) and the correction-to-scahng function /i(x) from Eqs. ( pT]) and ([l^) . Figure 
1 shows the scaling functions /o(x) and the correction-to-scaling functions fi{x) for n = 2 
and 3 in 3d. The amplitude of fi{x) is arbitrary. It is determined by the coefficient b 
introduced in ([T0|): the value 6 = 2 was used in Figure 1. The scaling functions and the 
correction-to-scaling functions do not show strong dependence on n and d for n > 2. For 
n = 1 and n > 2 the scaling functions are very different, especially in the small— x region. 
The reason for this is the presence of the sharp interfaces in n = 1 systems, which lead to a 
finite slope at the origin in /o(x) and a cubic small-x behaviour in /i(x) [Q. For n > 2, 
the small-x is quadratic for /q, and quartic for /i, with logarithmic corrections for even n. 

Within Mazenko theory the correction-to-scaling exponent u is non-trivial and depends 
on both n and d, with u < 4 for all n and d in nonconserved 0{n) models. The upper bound 
of 4 is obtained when d 00 (for any value of n) or n ^ 00 (for any value of d). 

A noteworthy feature of Figure 1 is that the correction-to-scaling function fi (x) is much 
larger than /o(x) at large x (the same feature was found for many of the models studied 
in [^). This means that, in fitting data, scaling violations at large-x should be given less 
weight in choosing fitting parameters (e.g. the scale length L{t)) than violations at small or 
intermediate x, because corrections to scaling are larger there. 



IV. THE ONE-DIMENSIONAL XY MODEL 



An exact solution of this model was first presented by Newman et al [22|. The solution 



yields an 'anomalous' growth law, L ~ t^/^, for the characteristic length exhibited by the pair 
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correlation function, compared with the usual L ~ t^/^ growth law of nonconserved models. 
Mazenko theory does not predict this growth law, for the simple reason that the theory 
has been built in such way that it might be expected to give qualitatively correct results 
only for systems with topological defects (i.e. n < d), since the n-component auxiliary field 
m(r, t) is defined in terms of the underlying defect structure. Despite this, the theory does 
a reasonable job of accounting for the behavior of systems with n > d, and in fact becomes 
exact in the limit n oo. However, systems with n = d + 1, which can support topological 
textures, are poorly treated by this approach. 

An exact solution for the nonconserved n = 2, d = 1 system is possible because the 
equation of motion for the order parameter becomes linear in the angle representation, 
= (cos ^, sin ^), which is natural in the hard-spin limit, where 0^ = 1. In this limit the 
free energy functional is simply F = (1/2) / dx{d(f)/dx)'^ = (1/2) / dx{d6/dx)'^. The zero- 
temperature equation of motion for model A, d(j)/dt = —SF/5(j), becomes, in the angle 
representation. 



dt dx'^ 



(17) 



which is a diffusion equation for the phase angle 9. Thus one characteristic length scale is 
the 'phase diffusion length', Lq = t^/^, but this is not the scale which characterises the pair 
correlation function. 

Equation (|T^) can be solved in Fourier space to give e^it) = ekiO)exp{-kH). In eval- 
uating quantities of interest such as correlation functions, one needs to specify the initial 
conditions. The probability distribution, P{[9k{0)]), for 9k{0) is conveniently chosen to be 
Gaussian P{[9k{0)]) oc exp — ^ /Sfc^fc (0)^-^(0) . The choice Pk = ^o^^/2 is made as it 
gives the initial condition C(r, 0) = exp(— r/^o) for the order-parameter correlation func- 
tion, which is the appropriate form for systems quenched from an equilibrium disordered 
state with correlation length ^q- The equal-time correlation function is given by 

C(r,t) = (0(x,t) ■0(a; + r,t)) = (cos[^(a; + r, t) -9{x,t)]). (18) 
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Using the Gaussian probability distribution for 9k{0) equation ([T8|), with the dynamics (|T 
gives 



C{r,t) = exp I -Y,-^exp{-2kH)[l- cos kr]] . 

\ k f^k J 



(19) 



Since the characteristic value of k in the integral is of order and we anticipate (see 

below) the growth law L(t) ~ t^^^, which sets the characteristic scale of r in C{r,t), it 
follows that the scaling limit and the corrections to it can be obtained from a power-series 
expansion of cos{kr) in (p!9|), since the characteristic value of kr is small (of order t~^^^) at 
late times. Retaining the leading and next-to- leading terms in the exponent, and evaluating 
the sums over k, gives 



C (r, t) = exp 



+ 



exp 



2eo(27rt)i/2 96^0 (27r) 1/2^3/2 

f2 „,4 
1 + ^ 



+ 
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L2 96(27r)i/2 ' - yL\ 

where y = r/L is the scaling variable and the coarsening scale L{t) = il'^t^'\ The 
correction-to-scaling exponent is = 2. 

This growth law is rather unusual since the generic form of the growth law for noncon- 
served fields is L{t) ~ t^/^. In this model u is found to be trivial while within Mazenko 
theory uo is non-trivial. There are two fundamental length scales in this problem, namely 
the phase coherence t^^"^ and the correlation length associated with the initial conditions. 
The coarsening scale L{t) of the pair correlation function is the geometric mean of these two 
lengths. Note that the pair correlation function has a strong dependence on which is 
not 'forgotten' at late times. This sensitivity to initial conditions is absent in other models, 
such as n = oo vector model, where the initial conditions drop out at late times. In the 
conserved 1-d XY model also, simulation results gives L ^ t^/s g instead of the L ~ t^/^ 



behavior expected in higher dimensions |T^, suggesting that this 'anomalous' behaviour may 
be present there also. Although no exact solution is known for the conserved case, heuristic 
arguments, based on the role of the two characteristic lengths, can account for the observed 
ti/6 growth ^|. 
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V. CONSERVED 0(N) MODEL 



The dynamical scaling properties of systems with a conserved order parameter (Model 
B) with 0(n) symmetry is studied using Mazenko theory. Naive application of this theory 
does not give correct growth law L ~ t^/^ for scalar fields (the bulk diffusion field must be 
included in order to get the correct law [^). Here, however, we will only consider systems 
with n > 2. For Model B systems, equation (HI) becomes 



1 dC{l2) 

2 dt 



V'C(12) + «(t)7 



dC(12) 
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(i7 

with a{t) = l/2S'o. For eq. (^) to have a scaling solution it is clear that a ~ 1/L^ and 
L ~ t^^^. The latter is the correct growth law for n > 2, but for n = 2, d > 2 there are 
logarithmic corrections which ( ^T]) fails to predict. We will first consider the case where 
n is very large. In this case an expansion in 1/n can be made in equation (^T]) . For large n, 
C(7) ^ ^ _ 7(1 _ ^^)/2n + 0{l/n^) and 7 = C + C^n + 0{l/n^). With the above 

truncations, fBll") can be written as 



l dC{12) _ 
2 dt 



V'C(12) + a(t) (c + — 
\ n 



(22) 



correct to order 1/n. 

It is worth mentioning that the C^/n term is essential for scaling to be recovered at finite 



n. For n strictly infinite the C^/n term is absent and 'multiscaling' is obtained [|Tl|]. For 
arbitrary n, an expansion in powers of C can be made. Truncating the expansion at order 

leads back to (0) but with n replaced by an effective n*, given by n* = {n + 2)a^ with 
an = n[5((n+l)/2,l/2)]V2 7r §. 

Dimensional analysis of (^2]) requires a{t) = a/L'^, which defines L. Including the leading 
corrections to scaling as usual we write 

C(r, t) = /o(r/L) + L--h{r/L) + ■■■ (23) 
dL/dt = + + ■ • • , (24) 
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where b fixes the amphtude of fi{r/L). Inserting these expansions into ( ]22|) and comparing 
terms of leading order, 0{1/L^), and next-to-leading order, 0(1/ L^'^'^'^^), gives 

Jo 



8 dx 

^ dx 8^^ 2 dx 



V^./o + « /o + 
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n 



(25) 
(26) 



where = d^/dx^ + [(rf — l)/x] rf/rfx. 

For general n one must solve (pi]) with C{r,t) given by (1). However, the singularities of 
(7(7) and its derivatives at 7 = 1 introduce some numerical difficulties. Instead, therefore, 
we solve ([2^) which is valid for large n. For general n, an expansion in C up to C^, leading 
to ([2^ ) with an effective n* 0], gives scaling functions which are in fairly good agreement 
with simulation results P,p!5[|. 

In solving (^5|) numerically, one must know the boundary conditions. These are provided 
by small-x and large-x analyses. For small-x, the series expansion /o = 1 + J^'i^iPrX^, 
substituted into (|2|), gives /o = 1 + /^x^ - (1 + 3/n) a(3x'^/A{2 + d) + with /^a = p. 
Numerical integration can therefore be performed on ( P^D with initial conditions /o(0) = 1, 
/g (0) = 2/3, /o(0) = /o"(0) = 0. Both a and P are undetermined parameters. 

For the large-x analysis, we impose the physical condition that fo{x) for x — >■ 00. 



This leads to the linearised version of eq. (^Sf) given by 



8 dx 



(27) 



There are four linearly independent solutions of (pT]), with the general asymptotic form 



fo{x) ~ Fqx'^ exp{—Bx'" — Ax'^ 



(28) 



The first solution is the constant solution, corresponding to A = c = B = 0. It satisfies (^7^ 
by inspection. The other three solutions are obtained by substituting (^) into (P7| ) and 
carrying out an asymptotic large-x analysis, leading to the relations 



v = A/3. 
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s = 2/3, 

= -l/8^;^ 
A = 64a5V9 , 
c = -2d/3 . 

The three different sohitions correspond to the three solutions for B, one real, two complex. 
The real solution, B = —l/2v = —3/8, leads to an exponentially diverging solution for /q: 

/o(x) ~ Fo x-'^'^l^ exp (3x^/78 - ax^'^), (29) 

while the two complex roots, B = 3(1 ± i\^)/16, generate two solutions which can be 
combined to give an exponentially decaying solution with oscillatory behaviour 

where Fq and (po are arbitrary constants. 

Just as in the Model A case, where A was fixed by imposing physical conditions on the 
large-x solution, also in this case we have an eigenvalue problem in which two parameters a 
and (3 are chosen to eliminate the unphysical constant solution and the exponentially diverg- 
ing solution. The same problem is encountered in Model B with a scalar order parameter 
P]. Applying the procedure described in it is possible to determine a and /?. 

Turning now to the corrections to scaling, we consider first the four linearly indepen- 
dent large-x solutions for the linearised form of eq. (|26|) . These are a power law solution, 
fi{x) ~ x~^, an exponentially growing solution, ~ x^exp(3a;^/'^/8 — ax^/^), and two decaying 
solutions that can be combined in the form 

3x^/3 / 373x^/3 ^VSx^/^ \ 
+ cos — + — + , (31) 

16 2 ; V 16 2 ^^ly ' ^ J 

where is arbitrary, andp = {uj — 2d)/3 ii u > A andp = {A — 2d)/3 otherwise. The small-x 
solution is /i (x) = /ix^ — a /i (1 + 3/^) x^/4((i + 2) + • ■ • . Therefore ( pB]) is solved numerically 
with initial conditions /"(O) = 2fi, /i(0) = /{(O) = /{" (0) = 0. The two parameters /i and 
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fo{x) 



-2d/3 



exp 



/i (x) ~ x^ exp 



uj are as yet undetermined. They are fixed in the same way as a and /3, by requiring that 
an oscillatory, exponentially decaying solution is recovered as x — cxo. 

Values for a, /5, and lo in 3-(i for n = 2, 5, 10, 20 and 50 are shown in Table III (For 
n = 2, the effective n* = n^/i has been used). The functions /o(x) and fi{x) are displayed 
in Figure 2 for n = 2 and 20 in 3-d. Again b has been set to 6 = 2 without loss of generality. 

The most important result to be extracted for Table III is that the value of u decreases as 
n increases. This behaviour is quite different from Model A, where u increases asymptotically 
to 4 as n increases. It seems from Table III that u probably tends to zero for n ^ oo, 
although an analytical determination of the correction to scaling for n large but finite, 
analogous to the treatment of the leading scaling function in [0 , has not yet been realized. 

Comparison of the scaling and correction-to-scaling functions displayed in Figure 2 rein- 
forces a point made in connection with the nonconserved systems, namely that the correction 
to scaling become large (relative to the scaling function itself) at large values of the scaling 
variable x. As noted before, this suggests that in carrying out scaling analyses of data, more 
attention should be paid to small and intermediate values of x, where corrections to scaling 
can be expected to be (relatively) smaller, than to large x. Indeed, for the nonconserved case 
(Figure 1) the correction to scaling has its maximum at a point where the scaling function 
is already quite small (around 0.1). 



VI. SUMMARY 

Corrections to scaling associated with a non-scaling initial condition have been studied 
in 0{n) models within the Gaussian closure scheme of Mazenko. We have calculated both 
the correction-to-scaling function, /i(x), and the associated correction-to-scaling exponent, 
UJ, for both nonconserved and conserved fields. In both cases Mazenko theory suggests that 
u is nontrivial, depending on the nature of the dynamics involved, the dimensionality, d, of 
the system and the number of order parameter components, n. For nonconserved fields the 
value of u tends to the limiting value 4 for n — oo with d fixed, and for d ^ oo with n 
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fixed. In the latter limit, the Mazenko theory reduces PJ2^ to the OJK theory |^ and its 



generalizations , believed to become exact as d ^ oo ||2^ . 

The 1-d XY model is anomalous in that it exhibits a different growth law from the 
standard one for nonconserved dynamics, and the correction-to-scaling exponent is simple 
(cu = 2). In this model quantities of interest, such as the correlation function C{r,t), retain 
'memory' of the initial conditions even in the scaling limit. 

In studying the conserved 0{n) model, an expansion in 1/n was used which is valid 
for large n. This approach was used to find the correction-to-scaling function and the 

exponent uj for n = 5, 20 and 50 in 3-d. For n = 2, an expansion in C up to was made. In 
the latter case, a comparison between leading-order scaling results and simulations shows 
very good agreement despite the wrong growth law (i.e without the logarithmic corrections 
predicted for n = 2 [[l^). In conserved systems u decreases as n increases, raising the 
question of whether u ^ or approaches some limiting value as n becomes very large. We 
have as yet been unable to find fi{x) and u analytically in the limit of large but finite n - 
this remains an interesting open question. 

The main lesson for the analysis of experimental and simulation data is that corrections 
to scaling can be expected to be relatively small at small and intermediate scaling variable 
X {=r/L), suggesting that this region be given more weight than large x in fitting (or 
collapsing) data. 
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TABLES 

TABLE I. Exponent A within Mazenko Theory for model A. 



n 


2 


3 


4 


5 


d = 1 


0.301 


0.378 


0. 11-1 


0.133 


d = 2 


0.829 


0.883 


0.912 


0.930 


d = 3 


1.382 


1.413 


1.432 


1.445 



TABLE n. Correction-to-scahng exponent lo within Mazenko Theory for model A. 



n 


2 


3 


4 


5 


d=l 


3.976 


3.982 


3.990 


3.993 


d = 2 


3.928 


3.946 


3.961 


3.970 


(i = 3 


3.930 


3.945 


3.958 


3.966 



TABLE IIL Values for the eigenvalues a, /?, fj, and lo within Mazenko Theory for model B. 



n 


a 






UJ 


2 


1.54880435 


-0.3336250 


0.227495 


2.4613967 


5 


1.72743447 


-0.3179775 


0.225025 


2.0667992 


20 


2.01748270 


-0.3292250 


0.214515 


1.1290901 


50 


2.179330049 


-0.3487125 


0.206515 


0.5029987 
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FIGURES 



Figure 1. The scaling function fo{x) and the correction-to-scahng function fi{x) for non- 
conserved order parameter. Continuous and broken hnes correspond to d — n — 3 and 
d — 3, n = 2 respectively. 

Figure 2. Same as in Figure 1, but for conserved order parameter. Continuous and broken 
hnes correspond to d = 3, n = 20 and d = 3, n — 2 respectively (note the different scales 
for the x-axes in the upper and lower plots). 
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